Transverse force on a quantized vortex in a superfluid 
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Abstract 

We have derived an exact expression for the total nondissipative transverse 
force acting on a quantized vortex moving in a uniform background. The 
derivation is valid for neutral boson or fermion superfluids, provided the order 
parameter is a complex scalar quantity. This force is determined by the one- 
particle density matrix far away from the vortex core, and is found to be 
the Magnus force proportional to the superfluid density. We conclude that 

contributions of the localized core states do not change this force. 
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For a vortex moving in a superfluid there is a force transverse to its velocity which 
is the counterpart of the Magnus force in classical fluid mechanics [0. This was used in 
the experiment by which Vinen established the quantization of circulation round a wire 
in 4 He. The obvious generalization of the form of the Magnus force to the situation in 
superfluid helium is to take the force per unit length to Fm = p s K x (vy — v s ), where 
K is the circulation, vy is the vortex velocity, and p s is the superfluid density and v s is 
the superfluid velocity in the neighborhood of the vortex. Even for the case of 4 He there 
has been some controversy about the form of this transverse force, despite the experimental 
measurements of Vinen || and of Whitmore and Zimmermann [f3| at temperatures where 
there is an appreciable difference between p and p s . This controversy is partly based on the 
observation by Iordanskii j3J that excitations such as phonons are asymmetrically scattered 
by a vortex, and this should lead to a transverse force proportional to p n and to the velocity 
v n — v y of the normal fluid component relative to the vortex. The magnitude of this term 
is not immediately obvious, nor is it obvious whether this contribution should be added to 
the Magnus force as we have written it here, or to the expression for Fm with the superfluid 
density replaced by the total fluid density. It has been convincingly argued || that the 
correct form of the transverse force is 

F = pK X (Vy - V s ) +p n K X (v„ - Vy) , (1) 

which can also be written as Fm plus a term proportional to v n — v s . The problem of 
determining the correct form of the force is discussed in detail by Sonin [[J. 

For a fermion superfluid there appear to be additional complications introduced by the 
spectrum of low-lying fermion excitations associated with the vortex core. Volovik |7],|8|] has 
identified an additional contribution proportional to v n — vy associated with these excita- 
tions, which, in some estimates, cancels most of the Magnus force. However, measurements 
of the circulation quantization in the B phase of superfluid 3 He [[]] show a Magnus force of 
the expected size. In superconductors there seems to be little unambiguous evidence about 
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the existence or magnitude of this transverse force. 

In earlier work we have shown the close connection between the Magnus force and a 
Berry phase [|Hj, but our arguments were dependent on a detailed model of the vortex 
structure. In this paper we consider a single, isolated vortex forced to move through an 
infinite superfluid with uniform velocity vy by some moving pinning potential. We are able 
to make an exact evaluation of the coefficient of vy in the force which has to be applied 
to the vortex in terms of the integral of the momentum density round a large loop that 
surrounds the vortex, and this is essentially the term proportional to vy in eq. ([!]). The first 
step in the argument involves a perturbative calculation of the effect of the motion of the 
vortex in terms of the instantaneous eigenstates of the Hamiltonian. In the case of a system 
with no inhomogeneous substrate this expression can be written in terms of a Berry phase 
[ 11 1 . This expression can be written in terms of the derivatives of one- and two-particle 
Dirac density matrices, or, alternatively, as the commutator of two components of the total 
momentum operator. This can be turned into an integral over a surface (or a loop for the 
two-dimensional case) surrounding the vortex and at a large distance from it. From this we 
can conclude that there are no corrections to the term proportional to vy in eq. ([I]). 

For simplicity we work with a two-dimensional system, but the argument can readily 
be applied to a three-dimensional system. We assume that the system is homogeneous and 
infinite. The isolated vortex in the superfluid can be pinned to a position r by applying 
a potential YU ^( r o — r i) which acts on all the particles in the superfluid. For a static 
problem this potential can be arbitrarily weak, but when the pinning potential is made to 
move through the superfluid it must have a strength sufficient that the vortex does not get 
detached from the pinning center, either by tunneling with the aid of the Magnus force, or 
by acquiring enough energy from the phonon system. The potential has to be strong enough 
to to break the degeneracy of this broken symmetry vortex state and allow us to use a 
perturbative treatment of the velocity. The situation is somewhat reminiscent of the theory 
of the Stark effect in atoms, where a perturbative treatment of a weak uniform electric field 
is possible, despite the fact that the electric field connects the bound levels of the electrons 
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in the atom to a continuum of ionized states of the same energy. The pinning potential 
may also be very strong, as it is in the Vinen experiment, where the circulation is pinned 
to a solid wire. Our results are quite independent of the details of the pinning potential, 
provided that it is strong enough to allow perturbative methods to be used, and does not 
break cylindrical symmetry. 

In terms of the instantaneous eigenvalues E a (t) and eigenstates \^ a (t)), for which we 
choose phases such that (^ a \^f a ) = 0, the time- dependent solutions of the Schrodinger 
equation can be written as 

|tf(t)) = a^e-'f^'W^ait)) + £ a^Oe-'/'W* 7 *!*^)) , (2) 

where, to first order in the velocity, a a (t) = 1, and 

a p (t) = - f dt'{^^ a )e- i ^ Ea ~ E ^ dt "l h 

= -f dt'(Vp\vv ■ V^a)e- i $ t ' {Ea ~ E e )dt " /h . (3) 

Here Vo denotes the partial derivative with respect to the position r of the pinning potential. 
This gives the expectation value of the force on the pinning potential as 

F = -£/«<*«|Voif|*a) 

a 

+ J2U^^oHe i I t ^ H)dt ' /h V a f dt'e- l f' {E <*- H ^ dt "l\ v • V + h.c.|^ a ) , (4) 

where H is the total Hamiltonian of the system, which depends on r only through the 
pinning potential, f a is the occupation probability of the state a, V a is the projection 
operator off the state a, and h.c. denotes the Hermitian conjugate term. The first term 
on the right of this equation includes any forces on the vortex proportional to the normal 
and superfluid velocities, and we have not attempted to evaluate these. We concentrate 
on the second term, which gives those forces proportional to the vortex velocity. There 
are energy- conserving dissipative forces included in this expression, but, as with the Kubo 
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formula for longitudinal conductivity, their evaluation requires a careful limiting process, 
and we have not done this here. We concentrate on the contributions from intermediate 
states whose energy is different from that of the initial state, and it is these that give rise 
to the transverse force perpendicular to the vortex velocity. 

For the homogeneous system that we are considering the eigenvalues are independent of 
time, so eq. @ can be written in the simpler form 



f = - E /«<*«! V #|*a> + 



ihV 

V H- ?—v v ■ V + h.c. 



(5) 



E a — H 

Since VqH is the commutator of the operator Vo with H, the commutator cancels the energy 
denominator, and the part linear in the vortex velocity can be written as 



F x n = -ihv v J2fa ^< 



dx 



dy 



dy 



dXn 



(6) 



where n is the unit vector normal to the plane. The Berry phase associated with a closed 
loop in ro can be written as the integral over the area enclosed by the loop of Fm/Tivv- 
Equation (P) corresponds to the familiar form for the Berry phase |fLl|| . Since we can choose 
the wave functions in such a way that the dependence on r is entirely through (r — r ), the 
partial derivatives with respect to r can be replaced by a sum over partial derivatives with 
respect to the particle coordinates r,-. Upon thermal average, this expression can now be 
written in terms of the Dirac density matrices as 



a i,i \ 



-ihh - / / d 2 r[V x V'p(r', r) 



ihh- J J d 2 n J J d 2 r 2 [2V l x v^r^.r^n,^)]^ 



(7) 



where p and V are the one- and two-particle Dirac density matrices for the system, and the 
sum over i and j denotes a sum over all the particles in the system. 

The integral over the two-particle density matrix T vanishes. Because of the symmetry 
between rj and the integrand in this integral can also be written as 



[(Vj + V[) x (V 2 - V' 2 )r(r' 1 , r' 2 ; r 1; r 2 ) - p(r[, r 1 )p(r' 2 , r 2 )] r=r , ; (8) 

the replacement of T by T — pp is also obvious, since the integral over r x of grad 1 p(r 1 ) 
vanishes. For any value of r 2 this can be integrated over ri to give a correlation function 
between distant points which should vanish if the range of integration is sufficiently large. 

A more formal way of getting the same result is to say that the first line of eq. ([?D 
gives the transverse force in terms of the expectation value of the commutator of the x and 
y components of the total momentum of the particles. This commutator is a one-particle 
operator given by 

dip^ dip dip^ dip \ 



[P x ,P y ] = J J dxdyl^ 



(9) 



dx dy dy dx J 

where the ip\ ip are creation and annihilation operators for fermions or bosons. Since it is a 
one-particle operator, its expectation value is given by the one-particle Dirac density matrix. 

The integrand of the first term in eq. (|7|) can be written as half the curl of (V — V)p(r', r) 
evaluated at r = r'. Now divide the integrals over r and into a sum over finite areas labeled 
with an index a; the first of these areas will be centered on the position of the vortex, and 
the others are all well away from the vortex. Stokes' theorem can be used to write the result 
as 

f/w = E / v rfr • t( v " v >( r/ ' r ) W • (io) 

a Jo z 

This result is exact, and there are no contributions to the integrals from the neighborhood of 
the vortex core, since we have chosen the boundaries of the regions of integration to be well 
away from the core; any contributions from the core states or from the properties of the core 
must be reflected in the density matrices well away from the core. In particular, because 
localized states inside the vortex core, such as occur for an s-wave superconductor ||12|| , do not 



influence the one-particle density matrix far away from the core, their contributions to the 
transverse force must be zero. It should be pointed out that this property is not transparent 
in eq. (|7|), where both detailed forms of density matrices and explicit integrations are needed 
to recover it. 
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For a neutral superfTuid the integrand in eq. ( fLOD is just the momentum density, equal, 
by definition, to p s u s + p n u n , where u s and u n are the local values of the superfluid and 
normal velocities. If the circulation of the normal fluid is zero, this term gives F = p s K x vy, 
the Magnus force Fm [PlRJTOf. From the Galilean invariance of the problem we can deduce 
that vy should be replaced by vy-v s , but we cannot determine the coefficient of the term 
proportional to v„ — v s without explicitly putting in different background velocities for the 
normal fluid and the superfluid components. 

At first sight it is surprising that we can obtain an exact result for this problem, since 
there are very few problems to which quantum many-body theory gives exact answers, but 
there is a similar argument that can be used in classical hydrodynamics. If fluid is streaming 
past a fixed vortex with velocity v, one can calculate the momentum balance on a very large 
cylinder centered on the vortex. The net pressure force per unit length on the interior of this 
cylinder is p n K x v/2, the momentum flow out from the cylinder is — p„K x v/2, and these 
must be balanced by a force — p n K x v acting on the vortex core. It is by an elaboration 
of this argument that Barenghi et al. |J obtain their expression for the contribution of the 
flow of the normal fluid past the vortex. 

Our conclusion is that for any infinite, homogeneous, neutral superfluid with a scalar 
order parameter the transverse force has the form p s K x vy, and those mechanisms which 
have been suggested as giving corrections must actually give alternative ways of looking at 
the same thing. Because our argument is, like the classical argument quoted in the previous 
paragraph, basically an argument that balances the forces on the vortex core with forces 
and momentum flow at large distances for the vortex, there is no obvious generalization that 
allows for the effect of the disordered background potential which is normally important for 
real superconductors. 
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